Finite- Temperature Phase Diagram of the Hubbard Model 
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The finite-temperature phase diagram of the Hubbard model in d — 3 is obtained from 
^ , renormalization-group analysis. It exhibits, around half filling, an antiferromagnetic phase and, 

between 30%-40% electron or hole doping from half filling, a new r phase in which the electron 
hopping strength t asymptotically becomes infinite under repeated rescalings. Next to the r phase, 
a first-order phase boundary with very narrow phase separation (less than 2% jump in electron 
density) occurs. At temperatures above the r phase, an incommensurate spin modulation phase is 
indicated. In d = 2, we find that the Hubbard model has no phase transition at finite temperature. 

^ : 

5^ , The Hubbard model is the bare-essentials realis- 

■ tic model of electronic conduction, yet essentially no 
, knowledge has existed even phenomenologically on its 
' most frontal macroscopic feature, namely its phase di- 

H \ agram at finite temperatures. In this research, we ob- 
tain a finite-temperature phase diagram for the Hub- 
bard model in spatial dimension d = 3, from an ap- 
proximate renormalization-group calculation with flows 
£j ' in a 10-dimensional Hamiltonian space. This rich phase 
'— diagram, in the variables of temperature, electron den- 

■ sity, and on-site repulsion, exhibits, around half-filling, 



an antiferromagnetic phase completely due to electron 
hopping. At 30-40% electron or hole doping from half 
filling, a new r phase occurs with distinctive conduction 
property. In the neighborhood of the t phase, a phase 
separation so narrow that the jump in electron density is 
less than 2% occurs. At temperatures above the r phase, 
an incommensurate frozen spin modulation phase is in- 
dicated. In d = 2, no phase separation or other phase 
transition occurs at finite temperature in the Hubbard 
model, in contrast pH to the closely related, but less 
realistic, t J model of electronic conduction. 

The Hubbard model is defined by the Hamiltonian 
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where c\ a and c J(T are the electron creation and annihi- 
lation operators with spin a =| or J, at site i of a cu- 
bic lattice, < ij > indicates summation over all nearest- 
neighbor pairs of sites, and 



f' 'in 



c\ a c i!7 and n { = ri lT + (2) 



are the electron number operators. The terms in the 
Hamiltonian of Eq.(l) are, respectively, the kinetic en- 
ergy term, the on-site repulsion (U > 0) term, and the 
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chemical potential term included in order to study the 
system over its entire density range from zero to two elec- 
trons per site. 

The renormalization-group transformation is formu- 
lated |||4|] by first considering a d = 1 system. An 
exact renormalization group transformation can be for- 
mally written, 

< U1U3U5... I n I viv 3 v 5 ... > 



E 



< U1W2W3W4M5W6... I e ° n I viw 2 v 3 W4V 5 w e ... > , (3) 



where Ui, W2,v 3 , etc. represent the single-site states. 
Primes indicate the renormalized system. The transfor- 
mation given in Eq.(3) conserves the partition function, 
Z = Z' , but cannot be implemented due to the non- 
commutativity of the operators in the Hamiltonian. An 
approximation is used: 

(even \ even 

~PH(i - l,i) - (3H(i, i + Uj] Tr wt exp (-{3H(i - 1, i) - pH(i, i + 1)) 

even / even \ 

= Jl exp (-(3'U'{i - 1, i + 1)) ~ exp ^ -/3'W(i -1,<+1) = exp(-/3'W) , (4) 

where 

-jSHfi, 3) = ~t (4c iCT + cj ffCj ,) (5) 

-(U/2d) (mmi + n ]\ n n) + (m/ 2 ^) E ( ni + n i) ■ 
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Thus, the approximation consists in neglecting the com- 
mutation relations beyond segments of three consecu- 
tive unrenormalizcd sites. This approximation is effected 
twice (~) in Eq.(4), in opposing directions, hopefully 
with compensatory effect. The crux of the calculation 
is extracted from the third step in Eq.(4), 

Tru)2 e-/3W(1.2)-/3«(2.3) =e -/3 , W'(1.3) . (6) 

When written in terms of three-site (on the left) and 
two-site (on the right) matrix elements, this equation 
amounts to contracting a 64 x 64 matrix into a 16 x 16 
matrix. This operation is facilitated by block diagonal- 
ization of the matrices, using the conservations of parti- 
cles, total spin magnitude, total spin z-component, and 
parity, so that the largest blocks are 4x4 and 2 x 2 for the 
unrenormalized and renormalized systems, respectively. 
Thus, a renormalized Hamiltonian —j3'7i' is extracted. 
The closed form of —f3'H' is more general than Eq.(l), 
namely 



0H= - ^2 [tohi- a hj^0 + t\ (m-rhj-e + hi^aUj-a) + t 2 n i - a n j - (7 ] [c\ a c ]a + c] 
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where the hole operator is = 1 — and the electron 
spin operator at site i is 



where 1? ' aa > is the vector of Pauli spin matrices. The four 
hopping terms in the flow Hamiltonian [Eq.(7)] corre- 
spond to one electron hopping with or without the oppo- 
site spin electron present at the initial and final sites (two 
of these processes are related by hermitivity and therefore 
have the same hopping strength t\) and to two electrons 
simultaneously hopping from one site to a neighboring 
site. These four processes can be called vacancy hopping 
(to), pair breaking (ii), pair hopping (£2), and vacancy- 
pair interchange (t x ). For 

t =t 1 =t 2 , t x = J = V 2 = V 3 = V 4 = , (9) 

the flow Hamiltonian [Eq.(7)] reduces to the Hubbard 
Hamiltonian [Eq.(l)]. Thus, Eqs.(9) are the initial con- 
ditions of our renormalization-group flows. However, 
in general, the hopping strengths renormalize differently 
and the new interactions are generated under rescaling, 
so that the renormalization-group flows are in the 10- 
dimensional, K = (to, t\, t 2 , t x , U, fi, J, V2, V 3 , V4), Hamil- 
tonian space. 

The transformation is implemented in d > 1 by 
using the Migdal-Kadanoff procedure, so that K — 
/ f)R(fK) where b = 2 is the length-rescaling fac- 
tor, the function it is the contraction process specified 
in the previous paragraph, and / is an arbitrary bond- 
moving factor, set to yield the correct transition tem- 
perature of the Ising model (/ = 1.2279 and 1.4024 in 
d = 3 and 2). This renormalization-group transformation 
yields known information about quantum systems, such 
as, in d — 1, the absence finite-temperature phase tran- 
sitions; in d = 2, a conventional phase transiton for the 
Ising model, a Kosterlitz-Thouless transition for the XY 
model HQ , no phase transition for the Heisenberg model; 
in d = 3, ferromagnetic and antiferromagnetic phase 
transitions for the Heisenberg model, the antiferromag- 
netic transition occurring at a 22% higher temperature 
than the ferromagnetic transition, a purely quantum me- 
chanical effect 0. The 10-dimensional renormalization- 
group flows also conserve the particle-hole symmetry, 
given the map: 

to = t 2 , ti = t\ , I2 = to , t x — t x , J — J , (10) 
fx = -fi + U + 2dV 3 - 2dVi , U = U + 4dV 3 - 2dVi, 
V 2 = V 2 -2V 3 + V 4 , V 3 = -V 3 + V4, V 4 = V 4 . 

The global analysis of the renormalization-group flows 
yields the phase diagram of the system. We have thus 
obtained the global phase diagram of Hubbard model, 
presented here in Figures 1-3, where first- and second- 
order phase boundaries are respectively shown by dottted 
and full curves. The particle-hole symmetry [Eq.(10)] 
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dictates that the Hubbard model [Eq.(l)] phase diagrams 
be symmetric about fi/U = 0.5, which is seen in all of 
our results. 

Figures 1 are for U/t = 20. Figure 1 (right panel) shows 
the full phase diagram in temperature versus chemical 
potential. Figures l(left, middle panels) show the de- 
tails in temperature versus electron density and chem- 
ical potential, respectively. It is seen that an anti- 
ferromagnetic phase occurs around half-filling, purely 
due to electron hopping, since the Hubbard Hamilto- 
nian [Eq.(l)] does not contain an explicit antiferromag- 
nctic coupling. In fact, we traced the occurrence of 
this antiferromagnetic phase to the non-zero value of 
the pair-breaking strength t\. The antiferromagnetic 
phase is unstable to at most 10% hole or electron dop- 
ing from half filling. Between 30 to 40% hole or elec- 
tron doping, a r phase occurs in which the vacancy hop- 
ping strength to or the pair hopping strength ti [see 
Eq.(7)], respectively, renormalizes to infinity under re- 
peated renormalization-group transformations. Thus, 
for hole doping, under repeated renormalization-group 
transformations, t — > oo, J/t Q = 2,V2/t = 3/2, fi/t = 
6,^o = 0, Z7 -> oo,U/U = 0,Vi/U = 0. Symmet- 
rically, for electron doping, the overbarred variables of 
Eq.(10) have this behavior. In all other regions of the 
phase diagram, all hopping strengths renormalize to zero 
under repeated renormalization-group transformations. 
Near the r phase, a first-order phase transition (dotted 
curves) occurs, seen as a single curve in Fig.l (center) 
in terms of electron chemical potential and opening up 
into a coexistence region in Fig.l (left) in terms of elec- 
tron density. The latter shows the distinctive feature of 
this first-order transition, namely that it involves a very 
narrow phase separation, e.g., a discontinuity in electron 
density of less than 2%. This is similar to what is seen 
experimentally in lanthanide compounds. Q At temper- 
atures above the r phase, a sequence of antiferromag- 
netic and disordered phases is seen, at many tempera- 
ture scales [Figs. 1,2 (center)]. We interpret this as the 
presence of an incommensurate spin modulation phase, 
with a temperature- and (less strongly, by the alignment 
of the sequencing) density-dependent periodicity. Our 
renormalization-group transformation, with a commen- 
surate rescaling factor and a built-in approximation, acts 
as a spurious substrate potential which, at small incom- 
mensuration, registers the incommensurate phase and, at 
large incommensuration, disorders it. The incommensu- 
rate phase that we thus deduce is indicated in Figs. 1,2 
(left). The features described above were also seen in the 
simpler, less realistic, tJ model. [||J] 

As U/t is decreased, the first-order phase boundary 
moves with respect to the r phase. It is seen that, for 
U/t = 4.44 (Figs. 2), it actually abuts the boundary of 
the r phase and, for U/t = 0.8 (Figs. 3), it is on the other 
side of the r phase. 

We have thus calculated a finite-temperature phase di- 
agram for the d = 3 Hubbard model that is rich in phase 
transition phenomena. We have also repeated the same 
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calculation for d = 2. We find that no phase separation 
(in contrast to the tJ model J|,f|]) or other phase transi- 
tion occurs at finite temperature for the Hubbard model 
[Eq.(l)] in d = 2. 
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Figure Captions 



Figs. 1: Calculated phase diagram of the d = 3 Hub- 
bard model for U/t = 20. First- and second-order phase 
boundaries are shown with dotted and full curves, respec- 
tively. Fig. 1 (right panel) shows the full phase diagram, 
which is symmetric about ji/U — 0.5. Antiferromag- 
netic [a] , disordered [D] , and r phases are seen. In the r 
phases, the hopping strength to or ti renormalizes to in- 
finity, for hole or electron doping respectively. Above the 
t phase, a sequence of antiferromagnetic and disordered 
phases [Fig. 1 (center panel)] is interpreted as an incom- 
mensurate spin modulation phase [Fig. 1 (left panel)]. As 
seen in (left), the first-order phase boundary has a very 
narrow coexistence region. 



Figs. 2: Calculated phase diagram of the d = 3 Hubbard 
model for U/t = AAA. 



Figs. 3. Calculated phase diagram of the d = 3 Hubbard 
model for U/t = 0.8. 
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Figure 3 




